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Abstract 

Obtaining reliable trace element reference values in tissues and fluids from inhabitants of the European Union [3] 
relies on the availability of standardized and harmonized protocols for the statistical treatment of the data on trace 
element levels in general European populations. In this context, cadmium was measured in the blood (BCd) of 514 
Italian inhabitants from the Lombardy region and the results statistically treated and presented according to a 
procedure which includes: simple descriptive statistics and graphical analysis such as stem and leaf and box-plot 
representations (average BCd levels were 0.62 n g/1; geometric mean, 0,51 /zg/l; median, 0-50 fig/1; mode, 0.30 
/r.g/1; 95th percentile, 1.48 p.g/1; 5th percentile, 0.20 jzg/1); p-p plot, Shapiro-Wilk and Lilliefors tests for normality 
(the distribution of the data is closer to the log-normal distribution and inconsistent with the hypothesis of 
normality); analysis of variance (BCd increases from 20 to about 60 years and then decreases; it is influenced by 
smoking but not by body mass and alcohol consumption and it is higher in men than in women); and step wise 
multiple regression analysis (BCd is influenced by the number of cigarettes/day and the total dose of exposure, 
cigarettes/day multiplied by smoking years). Tentative reference intervals for BCd based on the log transformation 
of the data are 0.14-1.82 jrg Cd/1 (whole population); 0.16-1.94 fig Cd/1 (male) and 0.13-1.66 fi g Cd/1 (female); 
0.24-2.68 /i g Cd/1 (smokers); and 0.14-1.27 fi g Cd/1 (non-smokers). 
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* Corresponding author. 

0048-9697/95/S09.50 © 1995 Elsevier Science BV. All rights reserved. 
SSDI 0048-9697(95)04474-2 


PM3001192979 


Source: https://www.industrydocuments.ucsf.edu/docs/txwk0001 


C Roggi el at. / Sci. Total Environ. 166 (19V5) 235-2T3 


236 

1. Introduction 

Knowledge of the concentrations of trace ele¬ 
ments in body fluids and tissues of persons with¬ 
out occupational exposure is the key parameter in 
establishing trace element reference values for 
evaluation of trace element intensity of exposure 
and uptake in humans [1]. Unfortunately, the 
existing information on trace element concentra¬ 
tions in the general population is generally of 
poor scientific quality and unsuitable to establish 
reference values to be used in regulatory decision 
making [2]. This is due to the complexity and the 
multidisciplinary character of the problem which 
must focus simultaneously on the composition of 
the reference population, prevention of sample 
contamination, documentation of quality control 
of the chemical measurements and proper statis¬ 
tical treatment [3], Not surprisingly, only a few 
research groups have managed to establish reli¬ 
able baseline data on trace elements in tissues of 
the general population. These observations are 
confirmed by the work carried out in the frame¬ 
work of the TRACY project [4J, an international 
project with the Eureka status which aims at 
establishing reference values in body fluids and 
tissues of non-occupationaily exposed subjects by 
quality grading of published articles on trace ele¬ 
ments in human tissues. Intolerable differences 
are found for trace element concentrations inves¬ 
tigated by different authors which reflect the in¬ 
adequacy and the inhomogeneity of the proce¬ 
dures used in selecting sample conditions, preana- 
lytical, analytical and statistical treatments [5]. In 
particular, the absence of established harmonized 
statistical treatments and presentation of the re¬ 
sults makes it very difficult to compare the data 
obtained in populations living in various coun¬ 
tries. This represents a serious limiting factor in 
establishing reference values for different popula¬ 
tions On a standardised basis. Thus, in designing 
new experimental procedures concerning the de¬ 
termination of trace elements in humans with a 
view to establishing reference values, such as in 
EURO TERVIHT (Trace Element Reference 
Values In Human Tissues) [6], a minimum stan¬ 
dardization of presentation of the data and statis¬ 
tical treatment is mandatory. 


The objective of this paper is related to the 
establishment of reliable cadmium reference val¬ 
ues in blood (BCd) of Italians, with particular 
emphasis on the adoption of statistical treatments 
and presentation of the data as a model to make 
possible the comparison of trace element refer¬ 
ence values derived for different countries [2]. 

2. Materials and methods 

2.1. Study population 

The population sampled consisted of 514 sub¬ 
jects. 231 males and 283 females, aged between 
20 and 79 years, a random subsample of 767 
individuals from the municipality of Rovescala 
(Pavia province). Northern Italy. They were en¬ 
rolled in a project for the study of the influence 
of alcohol consumption on human health 
(FATMA project, Consiglio Nazionale delle 
Ricerche). They lived in a hilly and viticultural 
area without industries, far from highways and 
with low traffic intensity. Twenty percent of the 
individuals analyzed were employed in industry 
and handcraft (mainly in manufacturing, food 
production and building); 28% in agriculture; 5% 
in commercial and other activities. Housewives 
and pensioners represented 47% of the total 
population group. All subjects had no history' of 
Cd exposure as derived from a questionnaire 
completed by specialists and which supplied de¬ 
tailed demographic, socio-economic and medical 
data. The population group was divided by strata, 
such as smokers and non-smokers (no consump¬ 
tion of cigarettes); non-drinkers (never consumed 
alcohol), and drinkers (> n g ethanol/day). The 
daily intake of ethanol was estimated for each 
individual from the consumption of wine, beer 
and superalcohol. 

2.2. Blood sampling and analysis 

Blood samples were collected with a commer¬ 
cial polystyrene syringe with a stainless steel nee¬ 
dle [3] into polycarbonate heparinized vials previ¬ 
ously washed with 1% HN0 3 subboiling and Mil- 
lipore water. Samples were stored at +5°C for no 
longer than 7 days before analysis as previously 
described [3]. 
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Cadmium concentration in the biood was as¬ 
sayed within 1 week of collection by graphite 
furnace atomic absorption spectroscopy (GF-AAS, 
Perkin Elmer, 5 100) after chelation and solvent 
extraction with sodium diethyldithiocarbonate- 
methylisobutylketone (Na-DDC-MIBK) and in¬ 
jection of 200 ju.1 of the organic fraction into a 
graphite furnace (detection limit, 0.03 fig Cd/1). 
The yield of extraction was checked on human 
blood spiked with l09 Cd 2+ radiotracer (0.1 jttCi 
109 Cd/ml of blood or 0.04 ng Cd/ml) or on blood 
of rats injected intraperitoneally with 0.1 fid of 
109 Cd iH \ Carrier-free 109 Cd radiotracer, as chlo¬ 
ride, was supplied by the Radiochemical Centre, 
Amersham, UK. 

Quality assurance of blood Cd was obtained 
from the analysis of certified standard reference 
material for Lead and Cadmium in Blood (Bureau 
Communitaire de Reference (BCR)); samples no. 
194,195, and 196 with a certified Cd content from 
0.5 to 12.4 fxg Cd/l. 

2.3. Statistical analysis 

The following statistical treatments were ap¬ 
plied: 

(i) Simple descriptive statistics (for each combi¬ 
nation of independent variables, the mean, stan¬ 
dard deviation (S.D.), geometric mean (G.M.), 
geometric standard deviations (G.S.D.), median, 


mode. 5th and 95th percentile were computed) ns 
well as graphical methods of analysis (stem and 
leaf, box-plot in order to describe the empirical 
distribution of the data). 

(ii) p-p plot Shapiro-Wilk and Lilliefors tests 
for normality in order to test whether the data 
follows a normal or log-normal distribution. 

(iii) Analysis of variance (ANOVA) in order to 
assess the significance of the differences in the 
group means. 

(iv) Step-wise multiple regression analysis in 
order to study the relationships between the log- 
transformed BCd measurements in males/ 
females and age, smoking habits, alcohol con¬ 
sumption and BMI (ratio body weight/height 2 ). 

The significance of the elaborated models were 
assessed by analyzing changes of the square of 
the coefficient of determination (R z ) as reported 
in the next sections. The SPSS package was used 
to analyze the data. 

3. Results 

3. J. BCd analysis 

Table 1 summarizes the results of the control 
of the preanalytical factors as well as the analyti¬ 
cal quality control of BCd analysis of blood in¬ 
cluding leaching tests on materials for sample 


Table 1 

PreanaJyticaJ factors and analytical quality control in the determination of blood cadmium in the Italian population 


Step 

Test 



Result 


Sampling (contamination) 

Leaching from vials 



0.03 M g/I 



a% HNO,) 





Storage (loss) 

119 Cd-spiked human blood 


0.08% 


Yield of extraction (before 

l0) Cd-spiked human (or rat) 


95.4% 


insturmental analysis) 

blood byNaDDC 





Standard Reference Material 

BCd (X ± S.D., ptg/l) 


CV (%) h 



(BCR, Pb and Cd in blood) 

Certified 

Determined 3 

Within 


Between 




series 


scries 

1 % 

12.4 +0.5 

12.2 +0.2 

— 


— 

m 

5.37 ± 0.24 

5.27 ±0.12 

4.9 


2.9 

194 

0.5 +0.1 

0.48 ± 0.04 

6.7 


4.3 


3 20 determinations. 

b 13 repeated measurements. 
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collection and loss to the walls of the containers 
during storage; the yield of Cd extraction by 
Na2-DDC-MIBK; the blank of the entire proce¬ 
dure, and precision and accuracy of the method. 

Contamination and loss of blood during sam¬ 
pling and analysis can be considered negligible 
for our purposes. The coefficients of variation 
(CV) determined in certified reference materials 
samples 194 and 195 to assess the precision were 
4.9% (within series) as well as 2.9% and 1.3% 
(between series) on 13 determinations. These re¬ 
sults and the mean deviation between the two 
series of 4% at the lowest level of Cd (accuracy) 
are considered satisfactory. The blank of the en¬ 
tire procedure was 0.01 jug Cd/1. 

3.2. Statistical distribution 

Table 2 shows the analytical results of the 


concentrations of Cd determined in the blood for 
all subjects, for age, sex and smoking habit. Fig. 1 
displays the corresponding box-plot representa¬ 
tion for all the data. The average level of BCd in 
all subjects was 0,62 fig Cd/1. Higher concentra¬ 
tions were observed in smokers (0.97 fig Cd/1) 
than non-smokers (0.56 p-g Cd/l); in men (0.69 
fig Cd/1) than in women (0.56 fig Cd/1); and in 
older rather than in young individuals. The values 
of mode, median and mean (mode < median < 
mean) indicate an asymmetry in the distribution 
of the data as confirmed by the box-plot in which 
the extreme high values in the data set influence 
the evaluation of the mean. This suggests that the 
data set does not fit a normal distribution. 

The probability plots of the raw data (Fig. 2A) 
and of the corresponding transformed log-normal 
data (Fig. 2B) and the Shapiro-Wilk test for nor- 


Table 2 

Blood cadmium concentration in inhabitants of the Lombardy region (Pavia province) in Northern Italy by sex, age, smokers and 
non-smokers 


Strata 


Blood cadmium ( fi g/1) 







n 

X 

<T 

G.M. 

O.S.D, 

Median 

Mode 

95th 

5th 

Age 

Whole 

514 

0.62 

0.48 

0.51 

1.89 

0,50 

0.30 

1.48 

0.20 


20-29 

65 

0.44 

0.38 

0.34 

2.04 

0.30 

0.30 

1.36 

0.10 


30-29 

52 

0.50 

0.34 

0.42 

1.86 

0.40 

0.30 

1.41 

0.16 


40-49 

64 

0.64 

0.43 

0.55 

1.70 

0.51 

0.30 

1.47 

0.20 


50-59 

99 

0.69 

0.47 

0.56 

3.63 

0.57 

0.40 

1.80 

0.20 


60-69 

155 

0.69 

0.62 

0.55 

1.62 

0.57 

0.50 

1.76 

0.20 


70-79 

79 

0.61 

0.34 

0.54 

1.66 

0.60 

0.70 

1.10 

0.21 

Sex 

Males 

231 

0.69 

0.58 

0.56 

1.86 

0.50 

0,30 

1.87 

0.20 


Females 

283 

0.56 

0.38 

0.46 

1.90 

0-50 

0.30 

1.30 

0.12 

Smokers 

Whole 

125 

0.97 

0.59 

0.81 

1.82 

0.80 

0.50 

2.20 

0.30 


20-29 

18 

0.75 

0.58 

0.57 

2.13 

0.65 

0.70 

2.06 

0.10 


30-39 

17 

0.83 

0.40 

0.74 

1.66 

0.70 

0.50 

1.44 

0.30 


40-49 

23 

0.87 

0.57 

0.76 

1.66 

0.70 

0.50 

2.70 

0.32 


50-59 

29 

1.14 

0.57 

1.00 

1.70 

1.00 

0.70 

2.25 

0.37 


60-69 

27 

1.18 

0.74 

0.98 

1.90 

0.90 

0.60 

2.86 

0.27 


70-79 

11 

0.79 

0.33 

0.72 

1.62 

0.79 

0.70 

1.30 

0.26 

Non-smokers 

Whole 

389 

0.51 

0.38 

0.42 

1.74 

0.45 

0.30 

1.00 

0.19 


20-29 

47 

0.32 

0.17 

0.27 

1.86 

0.30 

0.30 

0.75 

0.09 


30-39 

35 

0.35 

0.15 

0.32 

1.58 

0.30 

0.30 

0.75 

0.10 


40-49 

41 

0.51 

0.25 

0.46 

1.58 

0.46 

0.30 

0.99 

0.20 


50-59 

70 

1.50 

0.24 

0.44 

1.70 

0.50 

0.40 

0.96 

D.19 


60-69 

128 

0.58 

0.54 

0.49 

1.74 

0.50 

0.50 

1.17 

0.20 


70-79 

68 

0.58 

0.34 

0.51 

1.62 

0.55 

0.70 

1.10 

0.20 
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Fig. 1. Box-plot representation for the complete data set 
(Table 1). 


Fig. 2. Probability plot of the raw data (A) and the corre¬ 
sponding transformed log-normal data (B) (for details, see 
text). 


mality show that the log-normal distribution of 
the test (R = 0.98, P = 0.02) fits the data much 
better than the normal distribution (R — 0.74, P 
< 0.001). The Lilliefors test for normality (Fig. 3) 
confirmed this finding. 

Interestingly, the stem and leaf analysis of the 
transformed log-normal data has indicated a bi- 
modal distribution of the data set (results not 
shown) suggesting the presence of a factor which 
significantly influences the BCd in our population 
group. 

3.3. Influence of sex, age and smoking on BCd 

Box-plot representations of the transformed 
log-normal data (not shown here) reveal differ¬ 
ences in the data according to sex and age. Lower 
values were found in females compared to males. 
An increased trend of BCd was observed in sub¬ 
jects 20-59 years old while BCd decreased in 
subjects 70-79 years old. The analysis of variance 
on log-transformed data showed significant dif¬ 
ferences (P < 0.0001) between males and fe- 



Fig. 3. Lilliefors test for normality. 
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Fig. 4. BCd as a function of age. 


males as well as subjects with age < 30 years 
compared to the class of 40 years old (Fig. 4). 

When the raw data of Table 2 concerning BCd 
in all smokers and non-smokers had been trans¬ 
formed to log-normal form, a significant differ¬ 
ence (P < 0.001) was observed between the two 
groups, BCd being lower in non-smokers. 

3.4. Multiple regression analysis and analysis of 
variance 

The relationship of the variable BCd to other 
variables such as age, sex, number of cigarettes/ 
day, total number of cigarettes (smoking years 
multiplied by cigarettes/day), alcohol consump¬ 
tion and body weight investigated by multiple 
regression analysis and ANOVA analysis gave the 
following results (Table 3): 

(i) The effect of the number of cigarettes/day 


on BCd {R 2 = 28.9%, P < 0.01) is apparently 
(but not statistically significant) strengthened by 
smoking age (R 2 = 30.4, P < 0.001). However, 
when the total number of cigarettes smoked is 
considered, R 7 becomes 31.1% (P < 0.001). A 
control taking into account only the total number 
of cigarettes smoked during the smoking years 
does not change the situation (R 2 = 29%, P < 
0 . 001 ). 

(ii) No effect of sex, alcohol, and BMI on BCd 
has been observed while the age {R z = 35.2%, P 
< 0.001) influences BCd, because it influences 
the total number of cigarettes smoked. 

4. Discussion 

Standardization and harmonization of statisti¬ 
cal treatments are crucial aspects in making pos¬ 
sible the comparison of the data on trace ele¬ 
ments in tissues of the general population living 
in different countries and in establishing refer¬ 
ence values [3]. In this context, cadmium has been 
determined in the blood of 514 Italians living in 
the Lombardy region, Northern Italy, following 
the approach previously reported [3] which in¬ 
cludes; the composition of the reference popula¬ 
tion, defined in detail according to a previously 
developed questionnaire; the control of preanalyt- 
ical factors (Table 1) carried out by means of a 
109 Cd radiotracer which was needed to determine 
the yield of extraction of the very low amounts of 
Cd present in blood where it is bound to high and 
low molecular weight components [7]; and the 
analytical quality control, which was ensured by 


Table 3 

Multiple regression analysis between BCd measurements and some variables (such as smoking habits and BM1 K ) using the normal 
and log-normal models 


Variable 

Parameter 






Normal model 



Log model 



0 (.y itr’J 

r-value 

/"-value 

0<X 10'/ 

c-value 

/"-value 

Cigarettes/day 

3.686 

3.74 

0.0001 

6.6176 

9.226 

o.ooon 

Smoking age 

0.576 

1.08 

0.279 

0.6743 

1.671 

0.0954 

Total dose exposure" 

0.0961 

2.31 

0.021 

— 

— 

— 

BMI b 

0.847 

0.10 

0.917 

— 

— 

— 

Age 


— - 

— 

1.3140 

6.168 

0.00001 


“Cigarettes/day X smoking years. 
h Weight/heiglit 2 . 
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the use of existing certified standard reference 
material (Table 1). Concerning the statistical 
treatment of the raw experimental data on BCd, 
the first step was the presentation of the data in 
the form of simple descriptive statistics. This was 
carried out according to the recommendations of 
the TRACY project [4] (Table 2) in order to 
describe the central tendency, dispersion and 
shape of the distribution of the values. Then, the 
next question was to decide whether the data set 
followed a normal or log-normal distribution. Box 
and p-p plots (Fig. 2) as well as the Shapiro-Wilk 
and Lilliefors tests (Fig. 3) for normality show 
that the distribution of the data is closer to the 
log-normal distribution and inconsistent with the 
hypothesis of normality. Thus, the analytical data 
of Table 2 were transformed to log-normal data 
which were the basis for drawing tentative refer¬ 
ence values for BCd (Table 4) as well as carrying 
out further statistical treatments; 

(i) Stem and leaf analysis showed a bimodal 
distribution, which suggests the existence of a 
factor which significantly influences BCd in the 
studied population. 

(ii) Variance analysis, which has shown the 
variation of BCd as a function of age, sex and 
smoking in agreement with a recent publication 
on BCd in Italian population groups [11], BCd 
increased from 20 to about 60 years and then 
decreased (Fig. 4). The tendency for Cd accumu¬ 
lation in the body with age is known [8]. It should 
reflect the lack of a homeostatic mechanism for 
Cd and the induction by the absorbed Cd of a low 
molecular weight protein, metallothionein, which 
is stored mainly in renal cortex [9]. In addition. 


Tabic 4 

Tentative reference intervals for BCd in Italians living in the 
Lombardy region 3 


Group of subjects 

Tentative reference intervals (/xg Cd/!) 

All 

0.14-1.82 

Male 

0.16-1.94 

Female 

0.13-1.66 

Smokers 

0.24-2.68 

Non-smokers 

0.14-1.27 


"Calculated from log-transformed raw data for BCd (Table 
2) according to the formula C M./(C.S,D,) ! (lower limit! as 
well as O.M. X (G.S.D4 2 (upper limit) [13]. 


Table 5 

Stratification of the population sampled according to the 
variables sex and smoking habit 3 


Group of 
subjects 

All 


Male 


Female 

n 

% 

n 

% 

n 

% 

Total 

517 

100 

232 

44.8 

285 

55.2 

Smokers 

134 

100 

87 

64.9 

47 

35.1 

Non-smokers 

383 

100 

145 

37.8 

238 

62.2 


“Chi 2 =29.4, df=\. 


the influence of smoking on BCd reported in the 
literature in different populations [10] is also con¬ 
firmed by our data (Table 2) while alcohol con¬ 
sumption and body mass (BMI) had no effect on 
BCd. Higher levels of BCd in smokers than in 
non-smokers, in men and in higher age Italian 
population groups have been reported by Alessio 
et al. [11]. 

(iii) Multiple regression analysis and analysis of 
variance — interestingly, the influence of sex on 
BCd, shown by bivariate analysis, has not been 
shown by multiple regression analysis. This is 
related to the different distribution of males com¬ 
pared to females in the smoking group (Table 5) 
which suggests that sex could represent a con¬ 
founding factor in the relationship of BCd and 
smoking habit. Moreover, BCd is influenced by 
the number of cigarettes/day and the total dose 
of exposure (cigarettes/day multiplied by smok¬ 
ing years) (Table 3). 

The present study shows that the frequency 
distribution of BCd concentrations was approxi¬ 
mately log-normal and so the geometric mean 


Table 6 

Derived tentative BCd reference intervals in general popula¬ 
tion groups living in different Italian cities 


Place of 
residence 



BCd (fig/I) 

n 

G.M. 

G.S.D. 

Tentative reference 
intervals 3 

Verona*' 

162 

0.54 

2.02 

0.13-22 

Milan" 

199 

0.695 

2.10 

0.16-3.06 

Brescia" 

404 

0.475 

2.195 

0.099-229 

Pavia (this work) 

514 

0.51 

1.89 

0.14-1.82 


"All subjects (Table 2). 
h Froro ClZ], 
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Table 7 

Derived BCd tentative reference intervals for inhabitants of the Pavia province in Northern Italy compared to the general Italian 
case Iist b 


Population BCd (pg/l) 

Smokers Non-smokers 

rt G-M. GS.D. Tentative n G-M. G-S.D. Tentative 

reference reference 

intervals intervals 


Pavia (this work)* 125 0.81 1.82 0.24-2.68 389 0.42 1.74 0.14-1.27 

Italian case list 444 0.79 1.86 0.23-2.73 835 0.37 2.08 0.086-1.6 


a AII subjects (Table 2). 
b Frctn [12]. 

(G.M.) and geometric standard deviation (G.S.D.) 
have been used for establishing tentative refer¬ 
ence intervals (Table 4). The tentative reference 
interval derived here for all subjects (0.14-1.82 
/a g Cd/1) is narrower than corresponding values 
for people living in other Italian cities [12] in the 
Lombardy and Veneto regions (Table 6). How¬ 
ever, small differences are observed when sub¬ 
categories of smokers and non-smokers are con¬ 
sidered (Table 7) which suggests that in establish¬ 
ing reference values for BCd, it is mandatory to 
take into account variations in the smoking habit. 
In any case, the importance of the toxicological 
significance of these differences remains an open 
question [II]. 
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